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Motivated by the organization of crosslinked cytoskeletal biopolymers, we present a semimicro- 
scopic replica field theory for the formation of anisotropic random networks of semiflexible polymers. 
The networks are formed by introducing random permanent crosslinks which fix the orientations 
of the corresponding polymer segments to align with one another. Upon increasing the crosslink 
density, we obtain a continuous gelation transition from a fluid phase to a gel where a flnite frac- 
tion of the system gets localized at random positions. For sufficiently stiff polymers, this positional 
localization is accompanied by a continuous isotropic-to-nematic (IN) transition occuring at the 
same crosslink density. As the polymer stiffness decreases, the IN transition becomes first order, 
shifts to a higher crosslink density, and is preceeded by an orientational glass (statistically isotropic 
amorphous solid) where the average polymer orientations freeze in random directions. 

PACS numbers: 82.70.Gg, 64.70.Md, 61.43.Fs, 87.16.Ka 



Introduction. Semiflexible polymers are macro- 
molecules whose behavior is dominated by their bending 
stiffness. Recent years have seen an explosion of interest 
in them because some of the most important structural 
elements of the cytoskeleton (e.g. F-actin) and the 
extracellular matrix (e.g. collagen) fall in this category 
In vivo, these filaments appear in the form of 
networks or bundles. A first step towards understanding 
the behavior of such complex aggregates is to study in 
vitro solutions of many filaments whose interactions are 
controled by crosslinkers 0, B 0|- There is currently 
significant activity on the mechanics and elasticity of 
crosslinked bundles Q . Previous attempts to theoret- 
ically describe the formation of ordered structures from 
disordered solutions use a generalized Onsager approach 
Q or a Flory-Huggins theory In both approaches, 
the filaments are modeled as rigid rods and the role of 
thermal bending fluctuations (finite persistence length) 
is neglected. 

In this Letter, we consider randomly crosslinked net- 
works of wormlike chains (WLCs) , which are character- 
ized by two parameters: the total contour length, i, and 
the persistence length, Lp 0] . Permanent crosslinks con- 
nect randomly chosen pairs of WLCs, such that parallel 
alignment of the two chains participating in the crosslink 
is enforced (see Fig. [T]). Such crosslinks can be realized 
experimentally with short actin-bundling proteins such 
as fimbrin The parallel alignment acts like an effec- 
tive Maier-Saupe interaction, giving rise to nematic or- 
dering in the gel with an unexpected strong dependence 
on the stiffness of the chain. If the persistence length, 
Lp, of the WLC is sufhciently large, nematic ordering is 
observed right at the gel point. The degree of the order- 
ing transition within the gel fraction is discontinuous as 
in most nematic transitions. However, the orientational 
ordering is mediated by the crosslinks, so that only the 
fraction of chains in the gel exhibit nematic ordering. 
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FIG. 1: Crosslinker aligning two wormlike chains at arc length 
s and s'. 

Since the gel fraction goes to zero continuously at the 
gel point, the nematic order parameter is also contin- 
uous at the transition. For smaller persistence length, 
we find a true first order nematic transition inside the 
gel phase with the distance from the gel point increas- 
ing with L/Lp. Our model also predicts an orientational 
glass, which is characterized by frozen random orienta- 
tions of the localised chains in the gel. In this phase 
the rotational symmetry of the system is broken for any 
realization of the crosslinks but is statistically restored. 

Model. We consider a melt of N identical semiflexi- 
ble polymers modeled as wormlike chains (WLCs) with 
contour length L and persistent length Lp in three- 
dimensional space. The Hamiltonian of the system has 
one part related to the bending stiffness of the WLCs and 
another part, Uev, ensuring excluded- volume repulsion 
[HI: 

Here k — 2LpkBT denotes the bending stiffness of a WLC 
and ti{s) — dsri{s) is its tangent vector at arc length 
s (0 < s < L) with |ti(s)| = 1. We introduce M per- 
manent random crosslinks which constrain the system 
in such a way that they fix the positions of the corre- 
sponding segments to overlap and their orientations to 
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be parallel or antiparallel. The partition function for a 
specific configuration of crosslinks reads: 
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Z{Cm) = (n'5(r.e(se) -r,a4))'5(t.,(se) -met.aO)), 

where rrie = ±1 and the average, (...), is over all polymer 
conformations with the Boltzmann weight exp(— /37i). 
We treat the constraints imposed by the crosslinks 
as quenched disorder, and our goal is to calculate 
the disorder-averaged free energy, F — — fcBT[lnZ], 
where [...] denotes averaging over all crosslink con- 
formations which are determined by the number of 
crosslinks, M, their positions and their polarities. Cm — 
{ie, je', Se, s'^;me}^i- Wc assumc that a disorder con- 
figuration with M crosslinks follows the Deam-Edwards 
distribution 11211: 
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The physical content of this distribution is that polymer 
segments close to each other in the liquid (uncrosslinked) 
phase, irrespective of their relative orientation, have a 
high probability of getting linked. The parameter /x^ can 
be thought of as a fugacity controlling the mean number 
of crosslinks per WLC: [M]/N is of order fi^. 

Order Parameter and Free Energy. As the number of 
crosslinks is increased to about one per chain, we expect 
a gel transition to an amorphous solid state with a finite 
fraction of WLCs localised at random positions. This 
spatial localisation due to crosslinking implies for suffi- 
ciently stiff chains a restriction also on the orientation of 
the chains. In a phcnomcnological ansatz, we model the 
probability of a single monomer segment, s, to be found 
at position x with orientation u as 



< 5(x - r(s))(5(u - t(s)) >oc e 



(1) 



Here a denotes the preferred random position of the 
monomer, ^ is the localization length, and e is a unit 
vector with 77 specifying the degree of orientational order- 
ing. If the nematic aligning interaction induced by the 
crosslinks is strong enough, then a nematic phase arises. 
In particular for a (uniaxial) nematic gel ±e is the global 
axis of orientation. If the aligning interaction is weak, 
then an orientational glass is expected and e is equally 
likely to point in any direction. For a given number of 
crosslinks, the effective strength of the nematic interac- 
tion is controlled by the persistence length Lp. For large 
Lp/L ^ 1, one crosslink per chain is enough to approx- 
imately fix the orientation of the whole chain, whereas 
for Lp/L <C 1 many crosslinks are required to achieve 
nematic ordering. Hence we expect a phase diagram as 
shown in Fig. ([2|) with /i^ controlling the number of 
crosslinks and LjLp the polymer flexibility. This phase 
diagram is born out by the calculations, as we now sketch. 



In order to calculate the disorder averaged free en- 
ergy, we apply the standard replica trick, [In Z\ = 
lim„_+o ( [-^"l ^ !)/"-• We formally decouple individ- 
ual WLCs from one another by eliminating the single- 
polymer degrees of freedom in favor of collective fields 
[l3|. In the saddle-point approximation, the density of 
the system is uniform provided that the excluded- volume 
interaction is strong enough to prevent it from collapsing 
due to the crosslinks. Note that this interaction cannot 
induce an isotropic-nematic transition a la Onsager. The 
replica free energy per chain then reads: 



/{f^(k,u)} 



J2r / |f^(k,ii)p-l0gZ, (2) 



where 



z= (exp {^^^vY,^J^n{k,u) 



— ^ / dse-''-^(^)x5(u-mt(.)))) .(3) 

Here k = (kO, k^, k"), u = (u^, u^, ...u"), and 
denotes a sum over all wave-vectors except for k = 
(0, 0,...,0). The average is over an {n + l)-fold repli- 
cation of the single WLC Hamiltonian. 

In saddle-point approximation, the field satisfies the 
self-consistent equation: 
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(e-*-^'(^)j(u-mt,(s))) ,(4) 



where the average refers to the single chain "partition 
function" of Eq. ([3]). S7(k,u) acts as an order parame- 
ter which distinguishes between various phases, such as 
liquid, crystalline and amorphous solid with or without 
orientational order. As discussed above, we focus here on 
two scenaria for orientational ordering: nematic gels and 
statistically isotropic amorphous solids (SIAS), where the 
orientation of the chains is frozen in random directions. 
Orientational ordering is mediated by the crosslinks, af- 
fecting only chains in the same cluster. Chains in finite 
clusters, coexisting with the infinite cluster, move and 
reorient thermally. Hence only the fraction of localised 
chains, Q, i.e. those chains which are part of the in- 
finite cluster, exhibit orientational order. Therfore we 
generalise the order parameter of the isotropic gel to the 
following form: 



f)(k, u) (1 - Q) % Q + g uik, u) ^o,S2^ok° 



(5) 



The gel fraction is denoted by Q and macroscopic trans- 
lational invariance requires SJ^^pk" = 0. 

Statistically Isotropic Amorphous Solid. In the SIAS 
the orientation of the WLCs in the gel are frozen, such 
that the prefered direction fluctuates randomly from 
chain to chain. Rotational symmetry is broken for any 
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single realization of disorder but is statistically restored. 
(This type of order has also been predicted for other sys- 
tems da, [IBl)- The order parameter involves an average 
over all chains, which is equivalent to an average over all 
directions. Motivated by the phenomenological picture, 
we make the following variational ansatz for the order 
parameter 



Nematic Gel. The alignment of the crosslinked poly- 
mer segments may trigger an isotropic-to-nematic phase 
transition, provided a macroscopic cluster of crosslinked 
polymers exists. For a uniaxial nematic the two direc- 
tions ±e are equivalent. Hence we sum over these two 
directions ±e in the pheneomenological ansatz ([T|) and 
generalize the order parameter to include nematic order- 



w(k,u) 
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where the last factor ensures proper normalization of the 
order parameter, / du a;(k = 0, u) = 1. 

Using this ansatz, and taking the limit n — > followed 
by cxD, N — > oo, N/V — const., the free energy per 
polymer reads: 
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(7) 



where we have kept only the leading order terms in Q, 
and rj. B2, and -B3 depend on I = L/Lp through 
the lowest moments of the WLC conformational proba- 
bility distribution. 

Stationarity of the free energy with respect to Q, f , and 
r] yields a continuous gelation transition at /i^ = 1 char- 
acterized by a nonzero gel fraction Q ~ — 1), inde- 
pendent of /. These results are universal for the gelation 
transition in the saddle point approximation and have 
been confirmed for various models (e.g. [H, Q El)- 



The localization length and the degree of orientational 
order depend on the stiffness of the WLC. In the flexible 
limit, the localisation length is determined by the radius 
of gyration (^^"^ _ \')LLp and the orientational or- 

der goes to zero as rf' ~ Lp/L. In the stiff rod limit 
{I 0) ry2^2 (diverges at /i^ = 1. Our perturbative ap- 
proach breaks down in this singular case where a single 
increasingly long stiff rod is formed. 



w(k, u) = e 



sinhry 



cosh(7;u" -e). (8) 



The experimentally accessible nematic order parame- 



ter, 



N 



'5 = 4E^ / rf^(3(e-t,(.))2-l), (9) 
i=i -'^ 



is obtained from the generalised order parameter field as 



S = 



>ri f^^ , cosh (772;) o . 

' dx . . / / (Scos^a;-!) , (10) 



sinh(77) 



showing clearly that only the gel fraction contributes to 
the nematic ordering. Even if the nematic transition is 
of first order - as it will turn out - the jump in the ne- 
matic order parameter may be very small due to a small 
gel fraction and in fact S may even be continuous (see 
below). 

If we substitute the above ansatz, Eq. ([5]), into the 
saddle-point free energy, we obtain: 



f{e,v.Q)^ue,Q)+ur,,Q)+o{^) , (11) 



The positional part of the free energy, /p(f^,Q), is ex- 
actly the same as that obtained from Eq. (O by setting 
r] — 0, whereas the lowest-order in Q orientational part 
reads: 
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/n(.,Q) = ^ln{. 
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(cosh(?7 e • t(s)) cosh(77 e • t(s')))v 
sinh^ (77) 



^} ,(12) 



where (...)w denotes averaging over the WLC conforma- 
tions. For finite rj, stationarity of the free energy with 
respect to variations in Q and yields exactly the same 
results as in the gelation transition considered in the pre- 
vious paragraph. Since ^ ii^q term which couples 



positional and orientational localization is of higher or- 
der and can be neglected close to the gelation transition 

(g « 1). 

If we expand the nematic part of the free energy up 
to quartic order in 77^ , we obtain the typical Landau free 



FIG. 2; The first panel shows the schematic phase diagram. The other three panels show the free energy per WLC (in units 
of ksT) as a function of r] for I — 0.001 (second), I — 0.15 (third) and I — 0.4 (fourth) at fj? — fi^ — 1. 



energy which yields a first order phase transition from an 
isotropic phase with = to a nematic phase with finite 
77: 



a{lii^)T^^ + b{l,i?)i^^ + c{l,i?)Tf, (13) 



with h{l, /i^) < 0. The coefficients depend on the control 
parameters of the system which are the crosslink den- 
sity (through fi^) and the single WLC flexibility (through 
L/Lp). The free energy is shown for jj? = 1 and var- 
ious values of I in Fig. For stiff chains {I — 0.001), 
the global minimum occurs at finite 77 already at the gel 
point, whereas for more flexible chains {I = 0.15), the 
minimum at finite rj represents a metastable state with 
the global minimum at 77 = 0. In this case, the transi- 
tion to the nematic happens at a higher crosslink density: 
/i^ = 1 + 0{1). Finally, rather semiflexible chains (Z = 0.4 
in Fig. [2]) exhibit not even a metastable nematic state at 
/i^ = 1, but do so at higher croslink densities: a nematic 
metastable minimum appears at /if > 1, which becomes 



global at the transition point > 



The resulting 



phase diagram is shown schematically in Fig. ([2]). 

Since the degree of orientational order, 77, is finite at 
the transition, one cannot make quantitative predictions 
based on the Landau expansion of the free energy. How- 
ever, we can improve it by expanding (...)w in Eq. in 
Legendre polynomials. The truncation of the power se- 
ries expansion becomes less reliable as the WLCs become 
stiffer. We have checked the behavior of the free energy 
keeping terms up to 77^*, and the qualitative features of 
its dependence on the crosslink density are robust. As 
we increase the crosslink density beyond /i^, there is a 
value at which /n(?7) becomes unstable for large 77. At 
that point, our perturbative approach collapses, since it 
is based on the assumption of the positional-orientational 
decoupling close to the gelation transition which is valid 
for finite rj. Also, at the stiff rod limit {L/Lp 0), /n(77) 
becomes unstable at = 1. 

Conclusions. In this Letter, we have shown how the 
geometry of the crosslinks together with the stiffness of 
the constitutent chains control the orientational order 
of the random macromolecular network which is formed 



upon gelation. Whereas the gel formation is solely con- 
trolled by the number of crosslinks, the orientational or- 
der induced by localisation is sensitive to both the num- 
ber of crosslinks and the stiffness of the WLCs. Rather 
stiff chains tend to exhibit nematic order right at the gel 
point with a continuous isotropic to nematic transition. 
More flexible chains exist in a phase with frozen random 
orientation right at the gel point and show a first or- 
der nematic transition only at higher crosslink density. 
Future investigations may extend this work to consider 
different crosslinking geometries; one important example 
is polar ordering induced by crosslinks. It would also be 
of interest to study the possibility of combined nematic 
and SI AS orderings. 
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